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Abstract 

We compute the corrections from two-photon and 7-Z exchange in parity-violating elastic 
electron-proton scattering, used to extract the strange form factors of the proton. We use a 
hadronic formalism that successfully reconciled the earlier discrepancy in the proton's electron 
to magnetic form factor ratio, suitably extended to the weak sector. Implementing realistic elec- 
troweak form factors, we find effects of the order 2-3% at ^ 0.1 GeV^, which are largest 
at backward angles, and have a strong dependence at low . Two-boson contributions to 
the weak axial current are found to be enhanced at low and for forward angles. We provide 
corrections at kinematics relevant for recent and upcoming parity-violating experiments. 
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Two-photon exchange corrections have recently been found to play an unexpectedly im- 
portant role in elastic electron-proton scattering. Despite being of 0{a) smaller than the 
Born amplitudes, 27 exchange effects have been shown to have a strong angular dependence, 
which significantly influences the extraction of the electric form factor in Rosenbluth sepa- 
rations. Such corrections were found to resolve a major part of the discrepancy between the 
electric to magnetic proton form factor ratio, G^/G^j, determined via the Rosenbluth and 
polarization transfer methods (see Ref. l|] and references therein). 

Elastic ep scattering has also been used to probe the strangeness content of the proton, 
through measurements of the strange electric and magnetic form factors. This is achieved by 
scattering polarized electrons from unpolarized protons, and observing the parity-violating 
(PV) asymmetry Apy = (ct/? — o-l)/{o-ji + cr^), where is the cross section for a right- 

(left-) handed electron. The numerator in the asymmetry is sensitive to the interference of 
the vector and axial- vector currents, and hence violates parity 

In view of the large 27 contributions to electromagnetic form factors, it is natural to 
ask what effect the exchange of two bosons (7 or Z) may have on the PV asymmetries. In 
particular, since the extracted strange form factors appear to be rather small [sl, these two- 
boson exchange (TBE) contributions could affect the extraction significantly. In this paper 
we compute the relevant corrections to the PV asymmetry arising from the interference 
between single Z boson and 27 exchange amplitudes (which we denote by "^(77)"), and 



'erence amplitudes (denoted by "7(^7)"). 
l|] , which allows a more natural implemen- 



between the one-photon exchange and 7-Z inter: 
We use the hadronic formalism developed in Ref. 
tation of hadronic structure effects in radiative corrections at low four-momentum transfer 
squared Q^, where PV electron scattering experiment are typically performed. 

In their seminal work on electroweak radiative effects, Marciano & Sirlin {4] computed 
the '~f{Z'-f) contribution at = 0, both at the quark level and at the hadronic level using 
dipole form factors. The 2^(77) contribution was calculated in Ref. ^ within a generalized 
parton distribution formalism, applicable at a scale of several GeV^. More recently, Zhou 
et al. 61] computed the TBE effects on Apv within a hadronic basis using monopole form 
factors. 

In the present analysis, we account for the finite size of the proton by using realistic 
electromagnetic form factors in the loop graphs, determined self-consistently from a global 
analysis of elastic cross section and polarization transfer data [3] including explicit 27 ex- 
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change corrections. Furthermore, while an overall, factorized correction was applied to the 
PV asymmetry in Ref. 6|, here we compute explicitly the individual TBE corrections to the 
proton and neutron (or to the sin^ ^^y-dependent and independent) terms in Apy. 

In the Born approximation, the amplitude for the weak current mediated by Z exchange 
is given by: 

where sin^ 61^^ = 1 - M^/M^ is the Weinberg angle, with Mw (Mz) the W (Z) boson 
mass. For the corresponding electromagnetic one-photon exchange amplitude, Ai^, we refer 
to Ref. jy] for details. The weak leptonic current is given by a sum of vector and axial- 
vector terms, = Me(5'y7^ + 9Al^'~1^)'^e^ where Qy = —(1 — 4sin^ 9w) and g'^ = +1 are the 
vector and axial- vector couplings of the electron to the weak current. The matrix element 
of the weak nucleonic matrix is given by = um^'^un, where the current operator is 
parameterized by the weak form factors: 

r'z = rFf+'-^F,^ + r^,G^^, (2) 

where q is the four-momentum transfer and M is the nucleon mass. Assuming isospin sym- 
metry, the weak vector form factors -Fj^f a proton target are related to the electromagnetic 
form factors of the proton and neutron -F72'" (at tree level) by 

= (1 - 4 sin^ 9w)F^'2 " ^^^2 - i^i,2 , (3) 

where are the contributions from strange quarks. The small factor (1 — 4sin^^^) 
suppresses the overall contribution from the proton electromagnetic form factors. The weak 
axial- vector form factor of the proton is given by = —GY + G\, where is the strange 
quark contribution. 

In the standard model the parity-violating asymmetry Apy receives contributions from 
products of vector electron and axial-vector proton currents, and axial-vector electron and 
vector proton currents. It can be written as a sum of proton vector, strange, and axial- vector 
contributions: 

Apy = - (-^^) {Ay + A + Aa) , (4) 
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where 



Av = gAp 



Ak sin^ 6w) 



C"red 

{^r^lVr^s _|_ ^r^ivr^s \ 
0"rcd 



(5a) 



Aa = g'v vW+W^ , (5c) 

Crod 

where Gp = vra/ (-\/2 sin^ 9wM^) is the Fermi constant, and a^cd = ^{G'^eY + '^{^'mY the 
reduced unpolarized cross section, with e = (l+2(l+r) tan^(6'/2))~^ the photon polarization 
parameter and r = Q'^/AM'^. 

The parameters p and k, in Eqs. contain higher order radiative effects, such 
as vertex corrections, wave function renormahzation, vacuum polarization, and inelastic 
bremsstrahlung, which have been calculated previously and are well known. At tree level, 
p = K = 1. Beyond tree level, p and n also contain contributions from the interference of 
the Born and TBE diagrams, which we denote by Ap and An, respectively. The form factor 
G^/ implicitly contains higher order radiative corrections for the proton axial current, as 
well as the hadronic anapole contributions [2|, |3| . At tree level, and in the absence of the 
anapole term, G^^ G^ above. 

The contribution of the ^'(77) and 7(^7) TBE corrections to the PV cross section can 
be written: 

Aa™'' = 2 3? [M^^M*z + {M^z + Mz^^M*^] , (6) 

where yVi^^ {M.^z) is the two-photon (^l-Z) exchange amplitude. Since the asymmetry 
Apv is constructed as a ratio of the PV cross section to the unpolarized (parity conserving) 
cross section, one also needs to include corrections to the latter from the interference of one 
and two-photon exchange amplitudes (denoted by "7(77)")- These have been computed in 
Ref. [l| within the current framework. 

The sin^ 6w dependence of the TBE corrections can be obtained explicitly by calculating 
separately the proton and neutron contributions of Eq. ([3]) to the PV asymmetry. In so 
doing the sin^ ^vi/-dependent and independent parts can be evaluated and the Ap and An 



4 



corrections determined from the vector part of A-py in Eq. (l5a|) : 



4^ 4^ -I- 4" 
+ 

j^p,tree ^p,trec _|_ ^n,tree ' 



= -^f ;^ ' ^ , (7b) 



where Ay""^ is the TBE contribution to Ay from the proton (neutron), and Ay^'''^'^'^'^ is the 
corresponding tree-level asymmetry, with Aa'^^"'"'^ the electromagnetic two-photon exchange 
contribution to a^cd- The contributions to Ap arise therefore from the 'y{Z'~f) and ^(77) 
corrections, as well as from the electromagnetic corrections 7(77). On the other hand, Ak 
receives contributions only from the 'y{Z'~f) and ^'(77) corrections. 

The calculation of the TBE corrections proceeds along the same lines as that of the 27 

n 

amplitudes in Ref. with the replacement of the 'jNN vertex function by in Eq. ([2]). 
As in Ref. we parameterize the form factors as sums of three monopoles, but take 
the proton form factors from the more recent global fit of Ref. [31, and the neutron form 

n 

factors from Ref. [8| . Since the main purpose of the P V experiments is to extract strange 
contributions to form factors by comparing the measured asymmetry with the predicted 
zero-strangeness asymmetry, in all our numerical simulations we set the strange form factors 
to zero, F{*2 = = G\. For the axial- vector form factor we use the empirical dipole fit, 
Ga{Q'^) = G'a(0)/(1 + (5V^i)^ where Ga{0) = 1.267 is the axial- vector charge, and the 
mass parameter = 1 GeV. Varying by 20% does not affect the results significantly. 

In Fig. [U we show the calculated Ap and An corrections for = 0.01, 0.1 and 1 GeV^, 
relative to the soft-photon approximation (SPA) of Mo & Tsai j^. (All of the following 
results will be relative to the SPA.) In the SPA the Z exchange is factorized, which results 
in a cancellation of the model independent infrared contribution to the PV asymmetry. The 
total Ap is ~ 1-2% over most of the range of e, increasing at small e and small Q^. At 
low values the 2' (77) piece largely cancels with the 7(77), so that the total is saturated 
mostly by the 'y{Z'-f) contribution. At = 1 GeV^ the signs of the Z['-f'-f) and 7(77) 
corrections change, and the '^{Z'-f) contribution decreases in magnitude. 

The TBE corrections to k, are somewhat smaller in magnitude, ranging from ~ 0.2 — 0.3% 
at = 0.01 GeV^, where the Z^'-ff) contribution is dominant, to less than 0.1% at = 
1 GeV^, where there is large cancellation between these. As seen in Eq. fl7bp . there is no 
contribution to An from 7(77). 
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FIG. 1: TBE contributions Ap and Ak, relative to the SPA 9|], as a function of e for = 0.01, 0.1 
and 1 GeV^. The various contributions are: 7(77) (dotted), .^'(77) (dashed), 7(.^7) (dot-dashed), 
and the total (solid). 

The dependence of Ap and Ak is illustrated in Fig. [2] for several scattering angles. A 
rapid variation is evident for ^0.1 GeV^, especially at backward scattering angles. 
This is most pronounced in Ak and leads to a change in sign at ^ 0.05 GeV^ at 6* = 150°. 
Because of the strong dep endence, estimates at = 0, by Marciano and Sirlin ^ and 
more recently in Refs. lQ|,llI| in a somewhat different limit to the = point, are in general 
not sufficient to obtain a reliable correction for the actual experiments. Approximating the 
TBE corrections at non-zero by their = values may therefore lead to errors in the 
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FIG. 2: TBE corrections Ap and A/t as a function of for fixed scattering angle Q = 10° (dotted), 
e = 60° (dashed), 61 = 110° (dot-dashed), and 9 = 150° (sohd). 

TABLE I: Ratio 5 = / A^y of TBE contributions to the proton asymmetry relative to the 

tree-level asymmetry (in percent) at the and scattering angle 9 of selected past and future 
experiments. The results are compared for different input form factors (empirical, dipole and 
monopole). 



g2 (GeV2) 
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Ref. 










empirical dipole monopole 


0.1 


144.0° 


L12J 


1.62 


1.52 


1.72 


0.23 


35.31° 


[13j 


0.63 


0.58 


0.84 


0.477 


12.3° 


M 


0.16 


0.15 


0.24 


0.997 


20.9° 


[15j 


0.22 


0.23 


0.30 


0.109 


6.0° 


m 


0.20 


0.16 


0.32 


0.23 


110.0° 


m 


1.39 


1.33 


1.52 


0.03 


8.0° 


[18j 


0.58 


0.47 


0.86 



extracted form factors. 

The above behavior is qualitatively reproduced if one uses dipole form factors (for either 

the Dirac and Pauli, or electric and magnetic form factors, with a dipole mass of 1 GeV) or 

fl 

monopoles (for the electric and magnetic, with a monopole mass of 0.56 GeV |6(), instead 
of the empirical ones Quantitatively, however, there are significant differences between 
the empirical and monopole results at large e, with the monopole results for Ap (A/t) being 
~ 30% (60%) larger in magnitude at e ~ 0.9, with larger differences at larger Q^. This 
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reflects the sensitivity of tlie loop integrals to the large-momentum tails of form factors (and 
hence nucleon structure effects), even at low . 

The sensitivity to the input form factors is more clearly illustrated in Table I, where we 
list the values of the TBE contributions to the proton PV asymmetry relative to the tree- 
level asymmetry (in percent), (5 = Apy^/Apy , at the kinematics relevant to several past anc. 



future experiments 
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13 



14 



15 



16 



17 



18|. The results for the empirical form factors [?] 



are compared with those using dipole and monopole form factors in the loop integration. 

Generally the effects using the empirical form factors are ^ 0.5% for most of the forward 

angle experiments, increasing to ~ 1.5% at backward angles. The results with the dipole 
brm factors are similar, tending to be ~ 5-10% smaller. With the monopole form factors 
gI, however, the values of b are some 50% larger than with the empirical, which suggests 

insufficient suppression of contributions from large loop momenta. The TBE corrections to 

Apv may therefore be overestimated using monopole form factors. 

The impact of these differences on the strange form factors is difficult to gauge without 



performing a full reanalysis of the data, since in general 



different electroweak parameters 



and form factors are used in the various experiments 12, ll3|, ll4, ll5|, ll6|, ll7L _ 18j | . However, 
as an estimate of the possible effects we have determined following Ref. 6| the quantity 
defined there as as a measure of the induced difference between the strange asymmetry 
extracted using the different form factors. With the parameters of this analysis we find, for 
example, differences of the order of 15% between the empirical and monopole form factors 
for the HAPPEX kinematics 14 L around 20% for the GO datum 15| in Table I, and over 
30% for the PVA4 kinematics 13|]. Although these values should be treated as indicative 
only, they clearly point to the need for a more detailed reanalysis of the data including TBE 
effects and a careful treatment of the form factors in the loop integrations. 

As discussed above, the effective axial-vector form factor in the axial asymmetry 
Aa is defined to include the anapole form factor, and higher order radiative corrections. In 
extracting the strange form factors from data. Young et al. [sl fit the effective without 
decomposing it into its various contributions. Alternatively, one may extract the anapole 
form factor from the G^^ by correcting for the radiative effects. In Fig. [3] we show the 
correction 6a to the axial PV asymmetry, Aa Aa{1 + Sa), at several values. The 
correction is of order 1-2% for > 0.1 GeV^, but increases rapidly for decreasing Q^, 
especially at large e, where it reaches 10% at e ~ 0.95 at = 0.01 GeV^. This behavior 
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0.08 
0.06 

0.04 
0.02 






=0.01 GeV^ 




=0.1 GeV^ 




=1 GeV^ 




FIG. 3: TBE correction 6a to the hadronic axial part of the PV asymmetry, for = 0.01 (dashed), 
0.1 (sohd) and 1 GeV^ (dotted). 

may be related to the faster vanishing of the axial Born asymmetry compared with the TBE 
contribution at large e. 

The corrections calculated here are presented in a form that can be straightforwardly 
applied to the Apy data. The values for Ap and An can simply be added to the exist- 
ing radiative corrections contained in p and k, taking care to subtract any partial TBE 
contributions that have already been included ^. Because the TBE effects are largest at 
backward angles, they will be most relevant for the SAMPLE experiment at Bates 12| . 
and for the backward angle run of the GO experiment [l5| at Jefferson Lab. For the for- 
mer, we find {Ap,Ak){Q^ = 0.1 GeV^ ^ = 144°) = (1.950 x 10-2,7.998 x 10"^), while 

0.23 GeVl e 



for the latter {Ap,Ak){Q^ 
{Ap,Ak){Q^ = 0.62 GeV^^ = 110= 



= 110°) = (1.473 X 10-2,-1.342 x 10^^) and 
1.261 X 10-2, 1.103 X 10-^). In addition, even 



though it is at forward angles, the considerably smaller uncertainties expected in the Qweak 
experiment 18|] at Jefferson Lab will require a careful treatment of the radiative effects. 
In particular, we find {Ap,Ak){Q^ = 0.03 GeV^,^ = 8°) = (3.755 x 10-^,-2.655 x 10-^). 
A reanalysis of the entire data set on strange form factors incorporating these effects is 
currently in progress. 
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